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ABSTRACT 

New simulations for teaching quantitative biological 
techniques are now used at Michigan Technological University. 
Traditionallyr such techniques work within a particular system and 
have the student assume certain initial conditions and employ 
appropriate constants. The computer generates time dependent data 
which are plotted. The student then repeats the process with changed 
conditions, observing the changed effect. .In the new method, certain 
initial conditions are established randomly by the computer program 
and the student generates data to ascertain the unknown quantity, 
just as if he were measuring on a real system. The procedure is as 
follows: 1) The computer randomly generates a system parameter which 
the student must evaluate. . The value is stored unrevealed to the 
student. 2) The computer generates data based upon the unknown 
constant, student determined input values, and the prepared 
simulation program, 3) The student uses the data to determine the 
unknown. 4) Finally, the student compares his derived value with the 
stored value. This report also provides a detailed simulation 
example, a procedural flow chart, and a list of other potential 
training simulation. (PB) 



ERIC 



REPRINTED FROM! 



Proceedings of a 



CONFERENCE ON COMPUTERS 

IN THE 

UNDERGRADUATE CURRICULA 



The Claremont Colleges 
Claremont, California 

June 18,19, 20,1973 



THE INTRODUCTION OF BIOLOGICAL MENSURATION 
TECHNIQUES THROUGH SIMULATION 



JAMES D. SPAIN 



Sponsored by 
The Institute for Educational Computing 

of 

The Claremont Colleges 
with support from 
National Science Foundation 



THE INTRODOCTIOH OF blOLOGICAL MENSOBATION TECHNIQOES THROUGH SIMOLATIOM 



James D. Spain 
Michigan Technological University 
Houghton , Michigan 49931 
(906) 487-2029 



The use of simulation techniques for presentation of basic biological principles has 
expanded tremendously over the past several years. In most cases the student is provided 
with a prepared simulation involving a particular system, for example, predator-prey 
interaction. The student then assumes certain initial conditions, employs appropriate 
constants; and the computer generates time depeiident data. This data is either plotted out 
by the computer, or the student plots the data by hand. The student may then repeat the 
simulation vith changed conditions, and observe tiie effect. By working with the simulation, 
and by obtaining some understanding of the mathematics involved, the student gains a deeper 
insight into the basic principles and concepts of biology. 

At -the 1970 Conference on Computers in the Undergraduate Curricula, held at Iowa City, 
Iowa, I reported on a variety of biological simulation which could be ) mplemented on the 
Olivetti ffiicrocoiaputers E1]#[2]. Since that time, the list of available programs for this 
equipment has been extended to include simulations of both the deterministic and stochastic 
type in practically all fields of biological sciences [3]. 

During the past academic year, we have employed simulations in a slightly different 
manner at Michigan Tech with considerable success. Rather than having thp student select all 
the initial conditions for a givc^ simulation, certain initial conditions are established in 
a random manner by the computer program. The student them employs the program to generate 
data with the objective of accertaining the unknown quantity, just as he would if he were 
performing a measurement on a real system. The simulation becomes a teaching aid, for 
training students in the performance of some quantitative technique in biology. The 
advantage of this approach is that the student is provided with an individualized set of data 
which he must employ to evaluate one or more system parameters that are unique to his set of 
data. In terms of quantitative procedures, this training should be just as effective as if 
he were working with a real system. 

The basic technique is as follows: 

1. The computer employs a random number routine to generate the particular system 
parameter which the student is assigned to evaluate. The value is stored for use 
in the simulation without being revealed to the student. 

2. Next, the computer generates a set of data based on the unknown constant, st-udent 
determined input valuesr and the prepared simulation program. 

3. The student then employs the data to determine the unknown parameter, perhaps using 
some graphical procedure. 

4. Finally the student obtains a printout of the unknown parameter to compare with the 
value which he has derived from the simulation data. 

As an example of a training simulation of this type, consider the following which was 
developed for use in the ecology course at Michigan Tech, and is designed to introduce the 
'"student to the technique of estimating a population by the mark and recapture technique. The 
simulation has been programmed for the Olivetti p-'602 laicrocomputer. The computer flow chart 
of the simulation is presented in Figure 1. 



l?ark and Becapture Simulaticn 

The basic mark and recapture method of estimating populations of organisms involves 
taking a moderately large sample, marking them in seme manner, ' and returning them to the 
original group. At some subsequent time, additional samples are withdrawn and inspected to 
determine the fraction of marked individuals. The basic equation (Peterson) for estimating 
the population is as follows: 
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p = - • M 

"Where p is the population estimate, S is the nuaber sampled, r is the number of originally 
marked organisms vhich were recaptured, and a is the total Dumber o£ marked organisms in the 
population* 

This process is simulated using the Monte-Carlo technique by generating a random number 
to establish the true population size. This random number is then stored, in the computer 
without being printed out. The operator assumes a certain nusb<3r of originally marked 
organisms, which he enters i^to the computer. The sample size is also encored. The computer 
then sets up a ratio between marked organisms and the total population, dnd simulates the 
sampling process by generating a series of random numbers equal to the sample size. Each, one-: 
is checked against the ratic between marked and total. The computer keeps track of the 
number of times the random nufflber is greater than that ratio, and records these as unmarked 
organisms. The ratio is continually modified to reflect the changes in marked and unmarked 
organisms in the population as a ihole. At the end of the "sampling" the computer prints out 
the number of marked, and un narked organisms in the sample, the total number of marked 
organisms in the population, and :he population estimate based on three different equations: 

m (u-<-r ) 
r 



Peterson: p 
Schnabel: p 



Schumacher- 2 

Eschmeyer: p = gm (u+r) 

£inr 

where m = the total number of organisms marked. 

r = the number of marked organisms in the sample, 
u = the n'uraber of unmarked organisms in the saiple. 

At the end of each sampling the unmarked organisms are considered to have been returned to 
the population as marked organisms. 

The operator may have repeated "samplings" to observe the effect on the estimates which 
he obtains. When he is satisfied that his estimate is correct, he may call for a printout of 
the actual population. 

Figure 2 shows a plot of typical data obtained when 100 individuals were "marked" and 
then 11 samples of 50 individuals each were simulated by the computer. Note how the 
population estimates converge on the computer generated population ,of 1078. The stochastic 
nature of the simulation demonstrates quite nicely the magnitude of the deviations that one 
observes in this kind of measurement technique, and illustrates the peda-gogic value of the 
simulation* 

Although this approach has been employed for the development of only a few simulation of 
mensuration techniques, the following list suggests the broad applicability of the method in 
other areas of biological sciences. 

QtilSI EfitSBii^l Training Simulations in Biology InSiiJ^S- 
1. Delury estimate of population size. 
2m Michaelis-Menton constant for enzyme kinetics. 

3. Inhibition type, inhibition constant for enzymes. 

4. Titration curve for amino acids, pka determination. 

5. Extracellular fluid space determination. 

6. Molecular weight determination from diffusion data. 

7. Molecular weight determination from osmotic pressure. 
8« Half life and decay constant for radioactive decay. 
9« Chemical equilibrium constant determination. 

10. Chemical reaction rate constants, reaction orders. 
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11, Linkage strength and genetic mapping. 
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FIGURE 1, FLOWCHART OF THE BASIC 
MARK^ AND RECAPTURE SIMULATION 
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R « uniform random number 
P = population 
M = marked organisms 
s = sample size 
r = recaptured organisms 
U =» unmarked organisms 
p = population estimate 
S,E.= standard error 
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FIGURE 2 - GRAPH OF DATA 'FROM MARIC 
MD RECAPTURE SIMULATION 
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